Alternative method to find orbits in chaotic systems 
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We present here a new method which applies well ordered symbolic dynamics to find unstable 
periodic and non-periodic orbits in a chaotic system. The method is simple and efficient and has 
been successfully applied to a number of different systems such as the Henon map, disk billiards, 
stadium billiard, wedge billiard, diamagnetic Kepler problem, colinear Helium atom and systems 
with attracting potentials. The method seems to be better than earlier applied methods. 

PACS numbers: 05.45 


To analyze chaotic a system is it often necessary to find the periodic orbits in a system. This is especially important 
when applying semi-classical methods to investigate quantum mechanical systems. The periodic orbits (cycles) yields 
a skeleton for the dynamics [1,2]. We restrict ourself here to flows with two-dimensional Poincare maps and to maps 
of a two-dimensional plane. 

i n The problem of numerically finding the existing cycles is not trivial, and different methods have been applied. The 
O'Nproblem can be formulated as finding the zero points of complicated two-dimensional functions. This can be solved by 
^Vising the Newton-Raphson method, which converges very fast when the initial point is close to the correct solution. 

^One problem with this method is that good initial points are very hard to find for long periodic orbits. If one does 
'^not find a cycle it is difficult to know if the cycle does not exist or if one did not have a sufficiently good initial guess. 

For some systems there exist special methods. In a dispersing billiard system it is possible to find all orbits using 
Q^a “rubber band” method. This is simply done by choosing randomly one point on each boundary according to a 
description of which boundaries the particle should bounce off, and then vary these boundary points until the length 
^of the line connecting them reach a minimum. This is a “safe” method as it always converges to a solution which may 
^Ljbe the cycle or a non-physical solution in the case where the cycle does not exist. This method however converges 
^E^very slowly for long orbits and is therefore not ideal. A similar method has been used by Biham and Kvale for the 
[^.stadium billiard [3]. 

In the case of the Henon map there exists a method proposed by Biham and Wenzel [4] which converges to a cycle 
LOwith a given symbolic dynamics description. This method has a rather slow convergence and for some parameter 
^^values the method will not converge [5], 

More complicated methods have also been applied to find cycles in Hamiltonian systems [ 6 ]. 

The method proposed here works as follows: Let 5 and 7 be defined as the well ordered symbolic values for the past 
’’Oand future respectively (discussed below for specific examples). Assume we want to find the cycle with the symbolic 
QString S = S 7 S 2 • • • s n = ■ ■ ■ SiSo ■ ■ ■ s n SiS 2 ■ ■ ■ s n ■ ■ ■. We then make the following calculations 

• Find the symbolic values (<),..for the cycle S. 

• Calculate (5, 7 ) for a chosen starting point (x,y) in the Poincare plane. 

• Estimate the direction of the 7 - and d-axis in the Poincare plane close to the point (x,y). 

• Move the starting point in the Poincare plane by comparing (6, 7 ) for this point with the (S, 7 ) values for the 
cycle. 


Implemented as a binary search this procedure converges fast to the correct starting point (x,y) for the periodic orbit. 

The 7 -axis should be chosen parallel to the unstable manifolds and the d-axis should be parallel to the stable 
manifold at the point ( x,y ), but even rather large errors here only give a slow convergence and the method is very 
robust. It is usually enough to assume that a fixed direction of the 7 - and h-axis is valid in the whole Poincare plane. 
In the neighborhood of a homoclinic tangency this has to be done more carefully. If the system is a scattering system 
the method converges equally fast as for a bounded system using a proper definition of 7 and S for an escaping orbit. 
The binary search is done by using parallelograms in the Poincare plane defined by two corners and compare the S and 
7 values of the midpoint of the parallelogram with the correct values. When new untested points are chosen as corner 
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in a new and smaller parallelogram, the symbolic values of these points have to be tested since the manifolds are not 
exactly parallel to the 7- and 6-axis. This makes the method slightly slower than a straight forward one-dimensional 
binary search. 

If we want to find a non-periodic orbit with some particular symbolic description S, for example a homoclinic or 
a heteroclinic orbit, this is done by using (6,7) for this orbit. For a generic orbit which does not emerge from nor 
converge to a periodic orbit, 7 and 6 are not rational numbers. 

As a simple application of the method we describe how it is implemented for the Henon map [ 7 ]. The Henon map 
is defined as 


27+1 = 1 - ax'? + bxt-i (1) 

The Poincare plane is here simply the phase space plane (27, 27+1). For the parameter a sufficiently large this map 
yields a complete Smale horseshoe [8]. If |6| = 1 the map is area preserving. In the case of a Smale horseshoe the 
symbolic dynamics description for a starting point (*0,2:1) is given by the symbol string S = ■ ■ ■ s_os_iso • S1S0S3 • • • 
with s t . defined as 


JO if Xt < 0 

J 1 if xt > 0 


( 2 ) 


If the map is not a complete repellor we may use the partition defined by primary turning points as suggested 
by Grassberger and Kantz [ 9 ]. Well ordered symbols and the symbol plane (6,7) for this map was introduced by 
Cvitanovic and coworkers [ 10 ] to describe all admissible orbits by determining a pruning front in this plane. The 
definition of 6 depends on the sign of the parameter b because for b < 0 the map only folds the phase space once while 
for b > 0 the map both folds the phase space once and then reflects it yielding a different ordering of the unstable 
foliation. Both for b < 0 and for b > 0 is 7 defined as 


wt +1 


f m if s t+ 1 = 0 

\ 1 - m if s t + 1 = 1 


Wi = S 1 


7 ( 5 ) = 0.W1W0W3... = y wt/2 t . 

(=1 

The symbolic past 6 is defined for b < 0 as 


W t -1 = 




w t if s t - 1 = 0 


, iv 0 = s 0 


- w t if s t -i = 1 

OO 

5{S ) = O.wqW—iW— 2 ... = y^wi-t/2 1 . 


t: I 


and for b > 0 as 


Wt -1 = 


u 


- w t if S t -1 = 0 


w t if s t -1 = 1 ’ 

OO 

6 (S ) = O.W0W-1W-2 ... = ivi-t/ 2 1 . 

t=1 

For a period n orbit S7S0 • • • s„ we have 

7 = 2 2 "/ ( 2 2n — 1) 

and 

6= 2 2 "/(2 2 ” — 1). 


Wo = So 


( 3 ) 


( 4 ) 


( 5 ) 


The simplest assumption is to choose the 6-axis parallel to the 27-axis and the 7-axis parallel to the 27+1-axis in the 
whole (27,27+1) plane. In Fig.l we find that this yields approximately the 6-axis parallel to the stable manifolds and 
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the 7 -axis parallel to the unstable manifolds close to all crossing points between the stable and unstable manifolds. 
The parameter values are for this example a = 3 and b = 0.3 which gives a complete Smale horseshoe repellor. 

We choose a large rectangle defined by a lower-left corner (.T m i n ,*J nin ) and a upper-right corner (x max , x' ma ^). We 
then calculate the values 6 and 7 for the middle point ((;c m i n + a: max )/ 2 , (a:J llin + .rJ nax )/ 2 ) and choose the new rectangle 
to be 
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In case of the two last choices we do have to check the new corner points and move the lower-left point slightly down 
if the 7 value for this point is larger than 7 C and move it slightly to the left is the S value is larger than S c , and 
correspondingly for the upper-right corner. 

If a trajectory escapes we let the last symbol w t be 0 or 2 depending on the number of preceding symbols with 
s = 1 and let 7 and S be given by a finite sum. This gives the symbolic value for the folds at the border the escaping 
region. Implemented on a computer is it faster to just compare symbols w t than to calculate the values 7 and S. 

Instead of the binary search one could use a secant method which converges faster if one has a nice smooth function. 
The points S and 7 vary as a monotone devils staircase type of function with x and y (continous if the repellor is 
complete), and numerics yields that the simplest binary seach (bisection method) is faster for these fractal functions. 
Nummerics indicate that a secant method could converge faster if the devils staircase function only has short horizontal 
intervals, but not much is known of the convergence properties of such functions. 

A straight forward implementation of this algorithm for the repellor in Fig.l gives the starting points with an error 
less than 1CU 30 for all 4720 cycles up to length 15 in around 5 CPU-minutes on a normal work-station. We have 
also examined the method for the parameter values a = 1.4, b = 0.3 which yields a strange attractor, and for other 
parameter values. 

This algorithm is simple to program and we can use this method for many other systems. For the dispersing billiards 
consisting of 3-disks and 4-disks the algorithm for the symbolic values (S,j) is given in ref. [11]. Assuming the S- 
and 7 -axis are along lines ±mx in the Poincare plane with ( x , <f>) the posisiton and angle of a bounce the method 
converges to the correct orbit both for the complete repelling system, for the pruned repellor and for the overlapping 
disk bounded system. The number m depends on the parameters but can usually be chosen around 1.0. In this 
problem it gives no problems if there are pruned orbits, because all manifolds are sharply broken or ends in a point 
and there are no homodinic tangencies. 

The chaotic Wedge billiard [12,5] also yields a similar foliation of manifolds. The values of 7 and S are given by eq. 
(3) and (4) and the method works without problems. 

The Bunimovich stadium billiard has a more complicated symbolic dynamics. The well ordered symbolic values 7 
and S are defined in ref. [13]. Assuming 7 - and d-axis fixed at ±mx as for the disk systems, gives a working method 
to find the periodic orbits. 

A physically interesting system is the chaotic colinear Helium atom [14]. A Poincare map of this system using the 
point (ri,pi) each time ro = 0 (with ri,pi,r%,p 2 position and momentom of the two elektrons, see definitions in ref. 
[14]) yields a plane with manifolds ordered as a binary Cantor set and well ordered symbols are obtained from eq. (3) 
and (4). This Poincare plane is unbounded since there exist long cycles with arbitrary large r\ values, but this is not 
any problem when we only want to find orbits with a finite length. This system is a complete repelling binary Cantor 
set and we do not get any problems from homoclinic tangencies. 

The diamagnetic Kepler problem [15] (a classical Hydrogen atom in a homogeneous magnetic held) is a completely 
chaotic trinary repellor for sufficiently large scaled energies [16]. The orbits can be described by the same symbolic 
alphabet as the 4-disk system [17] and symbolic values are defined as for this system. For the complete repellor the 
symbolic dynamics is easy to determine and the method of finding orbits converges nicely using a similar Poincare map 
as for the colinear Helium. By numerically integrating this how we can hnd closing cycles in a few T CPU-secunds and 
starting points with 15 correct digits in a few CPU-minutes. We have calculated cycles as long as 1500 symbols, and 
we have applied this method in successfully calculating quantum mechanical resonances for Hydrogen in a magnetic 
held. For the pruned repellor case and for the bounded system one can describe the admissible orbits in this system 
using a pruning front description [18]. For these parameter values the symbols have to be dehned in a more careful 
■way and for some orbits one has to choose 7 - and d-axix with care. 

A different system where the presented method works is a Hamiltonian how in two dimensions with three attracting 
Gaussian potentials [19]. For different parameter values the dynamics can yield a binary Cantor set repellor or it 
may be a 6 -ary, 8 -ary, ... Cantor repellor and a symbolic dynamics can be given to the system. In the binary case we 
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define a “inner bent orbit” s = 1 and a “outer bent orbit” s = 0. Then eq. (3) and (4) yields the symbolic values and 
the method converges. A similar system with two Lenard-Jones potentials yields for some parameters a 5 Cantor set 
repellor and using a 5-letter well ordered alphabet we find periodic orbits in this system. The structure and symbolic 
dynamics of these systems will be presented in a later article. 

We expect this method to be applicable to all systems which can be described with a well ordered symbolic alphabet. 
If one wants to find an extremely long orbit, or an orbit with very high accuracy, then is it possible to find all points 
( x t ,yt ) with t £ {0,1,2 ..., (n — 1 )} by shifting the symbol string f-times and repeat the calculation finding ( x t ,yt ) 
by the method above. When integrating along such an orbit one adjusts the position at each crossing of the Poincare 
plane. We have however never found this necessary in any application. 

The author thanks for useful discussions Jan-Michael Rost and Predrag Cvitanovic, the leading cyclist. The author 
is grateful to the the Alexander von Humboldt foundation for financial support. 
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FIG. 1 
(5- and 7 - 



. Stable and unstable manifolds (TF S and W“)of the Henon map a= 3 and b = 0.3 and the estimated direction of the 
■axis. 
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